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As a model for mesoscopic quantum systems in thermal contact, we consider a four-mode Bose- 
Hubbard model with two greatly differing tunneling rates. By a series of Holstein-Primakoff trans- 
formations we show that the low-frequency dynamics of this system consists in general of two slow 
Josephson oscillations, rather than the single slow mode predicted by linear Bogoliubov theory. We 
identify the second slow Josephson oscillation as a heat exchange mode analogous to second sound. 



I. INTRODUCTION 

Experiments on mesoscopic quantum systems, and 
sophisticated numerical techniques for quantum many 
body theory, promise progress on deep questions about 
the relationship between microphysics and thermody- 
namics [IHB]. We still require analytical theory, however, 
to clarify the questions. For example, it should be well 
known that heat consists microscopically of energy held 
by degrees of freedom whose evolution is too rapid to per- 
ceive or control on macroscopic time scales [7] . This sim- 
ple formulation is nonetheless surprising to many physi- 
cists, because heat is rarely discussed explicitly except 
in terms of ensembles. Since the ensembles themselves 
are justified fundamentally as proxies for time averaging 
over rapid microscopic evolution, the role of time scale 
separation in defining heat is well established as an im- 
plicit assumption; but it is underdeveloped as an explicit 
theory. Investigating the mesoscopic onset of thermody- 
namics will require fully developing this theory, because 
the applicability of ensembles is not obvious in small sys- 
tems. 

Microscopic energy that remains adiabatically frozen 
in fast evolution may of course simply be ignored on 
longer time scales. What the First Law of Thermody- 
namics means from a microscopic point of view, however, 
is that energy can be exchanged slowly between slow and 
fast degrees of freedom: heat and work may transform 
into each other. To investigate the mesoscopic onset of 
thermodynamics, therefore, a simple dynamical model 
that provides a minimal representation of heat exchange 
offers a useful paradigm. 

In this work we present a simple but non-trivial model 
for two quantum subsystems in thermal contact, and 
show how heat exchange between them may be described 
in familiar quantum mechanical terms. Each subsys- 
tem is a pair of bosonic modes coupled so as to repre- 
sent an idealized Josephson junction - a two-mode Bose- 
Hubbard system. We specify that the subsystem char- 
acteristic timescale - the Josephson frequency - is to be 
the shortest one in the entire system, so that observa- 
tion and control of the system can address only much 
longer time scales. Each subsystem then has, in isola- 
tion, exactly two conserved quantum numbers: particle 
number and energy. For a fixed number of particles, the 
energy may still be varied independently by exciting the 



Josephson mode; if the rapid Josephson oscillations are 
averaged over, then states of any energy still count as 
equilibrium. Thus, the internal Josephson oscillations of 
each subsystem are to be regarded as heat. 

This identification is bold enough that it may need 
some defense. The 2-mode Bose-Hubbard model, which 
has been very extensively studied from many points of 
view, is certainly a rather minimalist model for a ther- 
mal system. It is technically ergodic, inasmuch as the 
corresponding classical system explores the entire energy 
surface (for fixed particle number) ; but only because that 
surface is one-dimensional, so that time evolution has no 
choice but to explore it all. The system is integrable, not 
chaotic; and chaos is often invoked as a justification for 
statistical mechanics. But the precise relevance of chaos 
to macroscopic behavior, as opposed to time scale sepa- 
ration alone, is far from fully understood. According to 
the explicit logic normally used to justify ensemble treat- 
ments of dynamical systems, and so relate thermodynam- 
ics to mechanics, the fact that time averaging reproduces 
the phase space average over the energy surface should 
allow us to consider the N-particle excited states of 2- 
mode Bosc-Hubbard as equilibrium states, as far as any 
measurements or control operations performed on much 
longer time scales are concerned. 

And if the standard appeal to time averaging allows 
us to claim our N-particle excited Josephson states as 
equilibria, then it is at least qualitatively reasonable to 
say that the energy of those excitations is heat. For a 
macroscopic sample of ordinary gas held in a fixed exter- 
nal potential, with a fixed number of atoms, there is a 
one-parameter family of equilibrium states, distinguished 
by total energy (usually expressed as temperature) . If we 
raise this energy without changing the external potential 
or the particle number, we say that we have added heat 
to the gas. Indeed this is almost the definition of heat, 
practically speaking. If heat means anything, it must 
mean this. 

We therefore argue that the excitations of a two-mode 
Bose-Hubbard system can indeed provide a simple model 
for heat, if they are involved in a dynamical process that 
occurs slowly compared to their Josephson frequency. We 
arrange for this by placing two such two-mode subsys- 
tems in thermal contact, through an additional Joseph- 
son tunneling between them, with a much smaller rate 
coefficient than that of intra-subsystem tunneling. The 
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two subsystems can thereby exchange particles, on the 
'slow' time scale; but they can also exchange heat, as 
the amplitudes of their fast Josephson oscillations beat 
slowly back and forth. Our combined four-mode Bosc- 
Hubbard system is thus a candidate minimalist model 
for heat exchange. 

The implications of this model for mesoscopic and 
quantum thermodynamics depend simply on whether our 
Josephson heat really behaves as heat. If it does not, be- 
cause our integrable system is in some way too trivial, 
then this will be a useful insight into the precise role of 
chaos. Even a reductio ad absurdum is still a reductio. 
But in fact we will show that the beating Josephson exci- 
tations do behave very similarly to heat in a recognized, 
though unusual, form: second sound, which is normally 
described as a temperature oscillation. We identify a 
nonlinear collective effect that makes Josephson oscilla- 
tions of heat - 'second Josephson oscillations' - the slow- 
est collective excitation of the 4-mode system about an 
excited equilibrium state. Our analytical derivation of 
this effect, from a simple quantum mechanical Hamilto- 
nian, indicates how thermodynamics may emerge from 
quantum dynamics in small systems. 



II. THE FOUR-MODE SYSTEM 



proximation, and show that the low-frequency behavior 
of the joint system consists of two Josephson oscillations: 
the elementary Bogoliubov excitation whereby the two 
subsystems exchange atoms with each other, and a non- 
linear beat mode by which they exchange heat. Insofar 
as a Josephson oscillation represents a two-mode analog 
of Bogoliubov zero sound in a dilute Bose gas, we identify 
our second Josephson mode as an analog of second sound 
(see for example [H]). We offer it as a simple example of 
a typical phenomenon of mesoscopic thermodynamics. 

Our system's four-mode Bose-Hubbard (BH) Hamilto- 
nian with onsite interaction U and index a = L,R is 



H = Hl + Hr — — (a\ L ai R + a\ L a 2 R + H. c 



H a — — — ( a\ a Cl2a + o-2 a a\ a 



(1) 



where ft is the high tunneling rate within each subsystem 
pair, and u) the low rate between the two subsystems. 

H conserves the total atom number N = Y^- aLa lr> , 
and in the present paper we assume for simplicity an 
eigenstate of N with large eigenvalue N. The global 
phase zero mode corresponding to N is therefore elim- 
inated, and we are left with the three nontrivial tunnel- 
ing modes illustrated in Fig. 2. In the linear Bogoliubov 



Our complete system is thus a rather symmetric four- 
mode Bose-Hubbard model, with two sharply differing 
tunneling rates. Such systems may be realized to a good 
approximation in current cold atom laboratories [9], 
either using four closely neighboring microtrap poten- 
tials, or with a double well trap holding atoms with two 
internal states coupled by a Rabi laser [THl [TT] . Our sys- 
tem can probably also be realized at least classically with 
superconducting Josephson junctions, but with the cold 
atom context in mind, we will refer to the particles as 
atoms. 
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FIG. 1: Bose-Hubbard system as model for thermal contact. 



By exploiting the small ratio between the slow and 
fast tunneling rates, we are able to analyze collective ex- 
citations in this model beyond the linear Bogoliubov ap- 
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FIG. 2: The three linear Bogoliubov modes. 

approximation the frequencies of these elementary exci- 
tations above the TV-atom ground state are 



L) = ^/uj(oj + UN), ft = v/ft(ft + UN) and 
ft' = ^/(fl + uj)(n + uj + UN). 



(2) 



These three modes are Josephson oscillations, and Q, ft, 
ft' are the three Josephson frequencies for the three in- 
dependent Josephson oscillations shown in Fig. 2. As is 
well known, these collective mode frequencies are shifted 
away from the pure tunneling (or for internal transitions, 
Rabi) frequencies w, ft by the nonlinear interactions. Re- 
pulsive interactions (U > 0) raise the Josephson frequen- 
cies, while attractive interactions (U < 0) lower them. 

In the oj <C ft case we consider, only one of these three 
elementary modes is slow, and the others are fast. Since 
excitations of the Bogoliubov modes provide a complete 
basis for the TV-particle many-body Hilbert space, it may 
seem puzzling that there are in fact (as we shall see) two 
distinct collective modes in this system with low frequen- 
cies (of order w rather than ft). This situation is actu- 
ally common in nonlinear quantum many-body systems, 
however. The 'type II' excitations in the one-dimensional 
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Lieb-Lieniger gas [131 HI > second sound in superfluid he- 
lium |12j . and even ordinary sound in air, are all examples 
of low-frequency collective excitations that are not to be 
found among a system's low-frequency elementary exci- 
tations, despite the fact that these are a complete set. 

The resolution of this paradox is that the frequency of 
beats between our system's two high frequency modes, 

ljj = (w/n)(fi + un/2) = o' - n + o(uj 2 /n) , (3) 

is of order uj, and it is the interference beat between the 
two nearly degenerate high-frequency modes that pro- 
vides the extra low-frequency collective mode. What we 
now show, by calculating to leading order in the small 
frequency ratio uj/fl, but nonlinearly in excitation ampli- 
tude, is that within the low-frequency regime this beat 
mode is an independent excitation in its own right, whose 
frequency is shifted away from wj by nonlinear effects. 



III. BEYOND BOGOLIUBOV 

To go accurately beyond the linear approximation 
within the standard symmetry-breaking Bogoliubov the- 
ory requires careful treatment of zero modes, but it is 
straightforward if we follow a number-conserving ap- 
proach [TS1 [TB] . Applying the Holstein-Primakoff trans- 
formation (HPT) [IT], we re-write H a= L :R m terms of 
the total atom number on each side N a = ^2 i=1 2 aL&ia, 
and atom-moving operators dj, , a a that commute with 
N ry : 



2\/ N a - ala a a a = (a\ a + a\ a )(a la - a 2a ) . (4) 

It is easy to show that [d Q ,dJJ = 1 and [a a , N a ] = 0. In 
these terms the Hamiltonians of the subsystems read 

H a = ~^N a + Qa\ a a 2a + 0(UN- 1 ) (5) 
+ jN a (N a - 2) + ^N a (al + a a ) 2 



U 



We can then implement the Bogoliubov formalism by 
defining a a = u a b a + v a w a and choosing u a ,v a to di- 
agonalize the quadratic terms in H a while maintaining 

[b a ,btt = 1. 

Since we are interested in dynamics slow compared to 
fl, and in Josephson oscillations of finite but not extreme 
amplitude, we will apply the rotating wave approxima- 
tion (RWA) by dropping from H a all terms that do not 
commute with the leading term proportional to b' a b a . As- 
suming the regime of large occupations for simplicity, we 
omit corrections of order UN^ 1 , and obtain for Hl r 



H a -> ~^-N a + ^N a (N a -2) 



(6) 



l m + 2UN a ) bib a - ^ n + 2UN « P3l. 

4 n + 2UN a 



In the quadratic term of ^ we note the familiar 
Josephson frequency of the excitations created by b\, 
which we have been unable to resist naming 'josons'. The 
co-efficient of the term quartic in b a is in general of the 
same order as U, but opposite sign; this reflects the well- 
known fact that the period of nonlinear Josephson oscil- 
lations lengthens with amplitude [19j . In quasiparticle 
terms, we may say that where atoms repel each other, 
josons attract; and vice versa. 

To simplify the description of the slow tunneling be- 
tween the L, R subsystems, we now perform a similar 
HPT for our two atom numbers, while assuming that our 
system is in an eigenstate of the conserved total Nl+Nr, 
with eigenvalue N. This defines a number-conserving op- 
erator d, which satisfies [a, a^] = I; it transfers atoms 
between the L and R subsystems. In the large N limit 
that we assume for simplicity, this provides 



N L . R = \[N±N 1 ' 2 {cJ +a)} 



(7) 



up to terms of C(d^d 2 iV _1 / 2 ). We can therefore re-write 
our full Hamiltonian |l| in terms of b^ r, a, and N. 

Since £1 is on the order of the problem's high frequency 
n, we see that to leading (zeroth) order in lo/Q,, our entire 
Hamiltonian is simply f2 J, for 

(8) 



J = b\b 



bW 



In classical terms, the 'total joson number' J is an adia- 
batic invariant. By means of an RWA we can now com- 
pute low frequency dynamics nonlinearly in excitation 
amplitude, but to leading order in small frequency ratios, 
by dropping all terms in H that do not commute with 
J. Since b]J) a are not separately adiabatic invariant, we 
must retain terms such as b^bp. Our result is 

ua)a + — N(a< + af 



H 



■^{b\b R + V R h) + ^ £ 

a=L,R 



i>ibib a b a 



Uj 



U flN 

TV n+uN 

n + UN/2 

y/n(n + un) 

U4fl + UN 

J n + uN ' 



(at 



a){b\b L 



(9) 



In the first two lines we can recognize two different forms 
of the standard single Josephson junction Hamiltonian. 
The first line provides Josephson oscillations of atoms 
betwee n L and R s ubsystems, with Josephson frequency 
uj = \Juj{uj + UN), while the second implies Josephson 
oscillations of josons. The two Josephson modes are ex- 
actly analogous, except that Uj and U are of opposite 
sign, as anticipated above following (5). The third line 
couples the atom and joson modes, and arises because 
the frequency of fast Josephson oscillations in each L, R 
subsystem depends on the respective atom number. 
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IV. SECOND JOSEPHSON OSCILLATIONS 

We can now perform a final HPT to express b a in terms 
of J and a J-conserving operator b that transfers josons 
between subsystems. As an adiabatic invariant, J is just 
as good a conserved quantity within the low-frequency 
regime as N itself, and so we may let J — > J in the low- 
frequency theory. To identify collective excitations, we 
then linearize (19| in a, b to obtain for fixed A, J 

H ViD = wtfa + ^N(a? + a) 2 

+ ojjPb + ^ J(6t + S) 2 

+ ~(NJ)^(tf + a)(tf + b). (10) 

In terms of the atom Josephson frequency ui — 
yJuj(uj + UN) and the new Josephson frequency for 
josons 



ojj = \/ujj(ujj + UjJ), 



(11) 



the decoupled collective mode frequencies in ( 10 ) are 

.1 l/ 2 



,~.2 



± 



~ 2 j\ 2 , ujujjW 2 NJ 



■ Ul 



n + uN 



(12) 



We observe that in the low-amplitude regime J <C A, 
either UN ^> uj so that uj ujj; or else UN < cj and 
thus uj fa uj j so that the third line in ( 10 ) is much smaller 
than the first two. Thus for J -C A the low frequency 
atom Josephson frequency (I> + = a; is always essentially 
unchanged from the standard Bogoliubov theory indi- 
cated in Fig. (2); but it is joined by a second Josephson 
oscillation whose frequency £>_ can be significantly lower 
(for U > 0), because Uj is negative. This is our main 
result. 



Note that the linearization that yields ( 10 ) is different 



from the standard Bogoliubov approximation of retain- 
ing only quadratic terms in fluctuations of the original 
a,i a about the mean-field ground state. It is a resummed 
Bogoliubov formalism, in which certain dominant nonlin- 
ear interactions among atoms have been included in the 
linearized interactions among josons, which are nonlin- 
early related to atoms through their definition via the 
HPT. One may also say that we do pursue the stan- 
dard (conserving) Bogoliubov approach, but applied to 
a two-mode system of interacting conserved atoms and 
josons. In this sense the extra low-frequency mode is a 
non- elementary collective mode in the full theory, but an 
elementary excitation within the adiabatic effective the- 
ory for low frequency dynamics. 

Although our derivation is post-Bogoliubov in the 
sense of being nonlinear, it is not essentially non-classical, 
because we have assumed large A, and because all our 
HPT and RWA steps are exactly paralleled in classical 
adiabatic theory. We can therefore test our theory of 



second Josephson oscillations numerically by solving the 
four-mode Gross-Pitaevskii nonlinear Schrodinger equa- 
tion for the classical limit of our system. We excite 
intra-subsystem Josephson excitations with an oscillat- 
ing drive, as follows. We relax to an A-particle mean 
field ground state, then excite josons by adding a time- 
dependent potential tilt to each subsystem, which varies 
periodically with the Josephson frequency f2. We ex- 
cite the second Josephson mode by giving the tilt drives 
slightly different amplitudes for the L and R pairs, so 
that slightly more josons are excited on one side than 
on the other. The resulting evolution, after the drive is 
turned off, is shown in Fig. 3 for a typical case. 




FIG. 3: Time evolution of AJ = b\b L - tf R b R and AA in the 
mean-field system for the parameters Q — 1, to = 0.1, UN = 
5, N = 10 4 and J = 139. Red solid line: numerical mean- 
field AJ. Blue solid line: oscillation at second Josephson 
frequency u~ computed in the text, fitted for phase. Blue 
dotted line: oscillation at beat-frequency ujj. Black dashed 
line: numerical AN. 

The agreement with our linearized theory is excellent, 
and the nonlinear shifting of the second Josephson fre- 
quency well below the trivial Bogoliubov beat is clear. 
Some atom population oscillation between L and R sub- 
systems occurs at frequency uj + , this mode having been 
slightly excited by our driving procedure, but the main 
atom oscillation is exactly opposite in phase to the jo- 
son oscillation. We note that while in incompressible 
supcrfluids second sound includes no include significant 
particle density change, the coupling of density and heat 
oscillations is expected for second sound in compressible 
superfhuds [20J. We interpret the coupling between the 
slow atom and joson oscillations as a simple example of 
the interconversion of heat and work. 



V. DISCUSSION 

The question whether our proposed oscillation of heat 
carries oscillating entropy is as much about the micro- 
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scopic definition of entropy, as about the validity of our 
identification of heat. Our josons do carry at least one 
physically meaningful form of entropy: the von Neumann 
entropy of the single-particle reduced density matrices 
Rl,r, projected onto cither half of our system. These 
are simply two- by- two matrices, diagonal in the basis of 
even and odd modes excited by a± a = (a la ± a2 a )/\/2: 



R a 



P a = 



a± a a±a) 
(N a ) 



I -Pa 

P a 



n + uN a 



+ 2UN a ) 



with 



1 \ 1 



(13) 



(14) 



for J a = (b a b a ). The local single-particle entropies S^r 
on each side, 

S a = -k B [P a lnP a + (1 - P Q )ln(l - P a )} (15) 

are genuine entropies inasmuch as they characterize the 
first order coherence, or equivalently the condensate de- 
pletion, within the two-mode Bose-Hubbard subsystems. 
They are direct analogs, for their respective two-mode 
subsystems, of the single-particle entropy which has been 
used since Boltzmann to represent thermodynamic en- 
tropy in statistical mechanics. 

Since S a is a simple function of J a , it is clear that our 
josons carry entropy at least in this single-particle sense. 
One could even consider letting J a — !• b\b a and N a —> N a 
in (14 1 and (151, in order to define Hermitian entropy 



operators S a which would be ordinary quantum mechan- 
ical observables. One of the basic questions in quantum 



thermodynamics is whether quantum entropy should be 
an observable operator just like quantum energy, with 
state-independent eigenvalues that appear as measure- 
ments with state-dependent probabilities, or whether it 
should remain a fundamentally statistical quantity like 
the full von Neumann entropy S = — fceTrpln p, which is 
a nonlinear function of the density operator p that repre- 
sents the quantum state of an open system. Our simple 
model thus addresses this question by suggesting a form 
of entropy operator, whose performance in the role of 
thermodynamic entropy may be assessed experimentally 
The general implication of our work for mesoscopic 
quantum thermodynamics would seem to be that it 
is rather easy to identify First Law phenomenology, 
whereby parametric oscillations among high-frequency 
internal degrees of freedom provide those extra contribu- 
tions to the low-frequency energy budget which in ther- 
modynamics are called heat. Inasmuch as this heat oscil- 
lates in our case, rather than diffusing, we must recognize 
that Second Law phenomenology, with irreversible spon- 
taneous processes, may be more elusive in small quantum 
systems. On the other hand, the full range of quantum 
dynamics in even the four-mode Bose-Hubbard model, 
in parameter ranges not considered here, is rich beyond 
the scope of this paper. It includes quantum chaos [2T] . 
whose role in mesoscopic irreversibility amply deserves 
further study. 
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